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The phase diagrams and transitions of nonequilibrium systems with multiplicative noise are 
studied theoretically. We show the existence of both strong and weak-coupling critical behavior, of 
two distinct active phases, and of a nonzero range of parameter values over which the susceptibility 
is infinite in any dimension. A scaling theory of the strong-coupling transition is constructed. 
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Though they have been argued jl| to describe a diverse 
and important set of physical systems out of equilibrium, 
among them autocatalytic chemical processes and several 
different problems in quantum optics, dissipative partial 
differential equations (pde's) with multiplicative noise re- 
main poorly understood. Such equations typically ad- 
mit two types of phases: trivial or "absorbing" phases 
in which the dynamical variable(s) or order parameter(s) 
vanish identically at all points in space, and remain zero 
in perpetuity; and nontrivial or " active" phases in which 
these variables have nonzero expectation values. To date, 
most of the theoretical effort on systems of this type has 
been devoted |l],|| to calculating the critical value of the 
control parameter at which the continuous transition be- 
tween these two phases occurs || . 

In this paper we present a somewhat broader investi- 
gation of the multiplicative noise problem, along the lines 
of existing analyses of critical phenomena in the directed 
percolation problem pj, and in the many physically- 
relevant stochastic pde's with additive noise. In particu- 
lar, we analyze the phase structure and critical properties 
of multiplicative-noise problems with different symme- 
tries. We show that above an upper critical dimension, 
d c — 2, the transition can belong to one of two distinct 
universality classes described by different fixed points: a 
weak coupling, mean-field fixed point, accessible for noise 
strength D less than a (nonuniversal) critical value D c ; 
and a strong coupling fixed point with nontrivial expo- 
nents, accessible for D > D c . The separatrix at D = D c 
is described by a third, unstable fixed point, likewise with 
nontrivial exponents. For d < 2, only the strong coupling 
critical behavior occurs. We analyze the mean-field tran- 
sition, compute the critical exponents on the separatrix 
in an expansion in e = d — 2, and formulate a scaling 
description of the strong coupling transition, which we 
check against previous exact results for d = (the single- 
variable problem). We also demonstrate the occurrence 
of a rather striking phenomenon: For discrete symmetry 
and any d, there is a region of the phase diagram in which 
the response of the system to an infinitesimal uniform 
field (i.e., the uniform susceptibility), is everywhere in- 
finite. For systems with N-component order parameters 
and O(N) symmetry with N > 2, we argue for the exis- 



tence of two distinct active phases for d > 2, one which 
breaks the O(N) symmetry and one which does not. We 
illustrate this with an exact calculation for N = oo. 
The models we study arc defined by the equation 



dtrii(x, t) = fjtV z Tii — mi — w|n| 2 / 'N + \n\rj. 



(1) 



Here it = {ni, ...,nzv} is an N-component, real vector 
field, r and u (> 0) are real parameters, and if is a Gaus- 
sian noise vector with correlations < rji(x, t)rjj(x', if) >= 
D8ij5(x — x')5(t — t')/N, for some noise strength D. We 
deal primarily with the cubic nonlinearity p — 1, though 
for the single-component theory, N = 1, we also discuss 
aspects of the quadratic nonlinearity, p = 0. We work 
exclusively in the Ito || representation. 

To establish some of the general phenomenology of 
these models, we first discuss mean-field theory. The 
neglect of fluctuations which is the essence of mean field 
is most simply implemented by dropping the noise term 
in (|l]), and ignoring spatial variations in it(x,t), which 
is replaced by it(t). For N = 1, Eq. (Q) then has two 
steady-state solutions with dn/dt = 0, respectively rep- 
resenting the absorbing and active phases: n(t) — 0, sta- 
ble for r > 0, and n = (—r/u) 1 ^ 1+p \ stable for r < 0. 
For the continuous symmetry case, N > 2, (where only 
the cubic nonlinearity, p = 1, preserves the O(N) sym- 
metry and so need be considered), a result identical to 
that for N = 1, p = 1 holds, except that in the active 
phase the orientation of n is unspecified, the magnitude 
being given by n 2 = {—r/u). In both cases the transition 
between the phases occurs at r — r c — 0, and the critical 
exponent (3 characterizing the vanishing of n as r — *• r c 
takes the value (3 = 1/(1 + p), independent of dimension, 
d. The correlation length and dynamical exponents as- 
sume their usual mean-field values, v = 1/2 and 2 = 2, 
respectively. 

Alternatively, one can formulate mean-field theory in 
terms of a model with infinite-range interactions (or 
equivalently, d — oo), by replacing p\7 2 n in (|l|) by 
a{< n > —it), where a is a positive constant and < n > 
the expectation value of n. For N = 1, this model can 
be solved self-consistently for < n > , giving a continuous 
transition at r c = with an exponent f3 = 1/(1 + p). 
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It is particularly easy to determine the critical dimen- 
sion d c above which this mean-field result is valid, by 
considering the absorbing phase where < n(x,t) >= 0. 
In this phase, one can calculate certain correlation and 
response functions exactly, by summing perturbation ex- 
pansions to all orders ||. We first specialize to N = 1, 
where the exact renormalized noise, Dr, and coupling 
constant ur, are given by Dr = DA, and ur = uA, 
with A = 1/(1 - DI d (r)), and I d (r) = J <0 533^7. 
Note that I d {r) is finite for all r > 0, provided d > 2 
. Hence for small enough D the denominator of A re- 
mains bounded below by a positive number, even when 
r decreases to zero. Thus Dr and ur remain finite and 
nonzero throughout the absorbing phase, right down to 
the critical point, which continues to occur precisely at 
r = 0, provided D is small enough and d > 2. In this 
regime, all other quantities of interest likewise experience 
no singular renormalization near r = 0, so the mean-field 
critical point at r = and mean-field exponents remain 
exact. (In fact some quantities, such as the response func- 
tion g (k, u>) — —iuj + pk 2 + r, experience no diagrammatic 
renormalization whatsoever in the absorbing phase ||.) 

For d < 2, by contrast, I d (r) diverges at r — 0, so Dr 
and ur both blow up for some positive value of r even for 
infinitesimal D. This invalidates mean-field theory, and 
suggests that the critical value r c is shifted away from 
zero, and that the mean-field critical exponents can no 
longer be trusted. The same is true for d > 2, provided D 
exceeds the special value D c for which Dr first diverges 
at r = 0. These results strongly suggest that d c = 2, 
but that even for d > d c there may be nontrivial critical 
behavior at large D 

This expectation is confirmed by straightforward 
renormalization-group (RG) analysis § of models (@), 
in the absorbing phase and at the critical point. Imagine 
rescaling space, time, and the field according to x = bx' , 
t = b z t', and n(x,t) — b^ri' '(x',t'), where z and C are 
as yet undetermined exponents. Writing b = e , and us- 
ing standard methods, and the results quoted above, one 
readily derives the following recursion relations 

dp/dl = (z — 2)p 
dr/dl = 2r 

du/dl = u((l + p)( + z + (1 + 2p)A d D) + 0(pu 2 ) 
dD/dl = D{z-2-e + A d D) . (2) 

Here e = d — 2, and A d — l/Aw + 0(e) is a posi- 
tive, d-dependent constant. Owing to the aforementioned 
absence of any diagrammatic renormalization of the re- 
sponse function in the absorbing phase, the p and r equa- 
tions are exact. The u and D equations contain no fur- 
ther diagrammatic corrections in the case p = 0. For 
p = 1, there are, as indicated in (^), corrections to the u 
equation of 0(u 2 ). 

Obviously, finite fixed points of these recursion rela- 
tions can only be reached by choosing z = 2 and tun- 



ing r(l — 0), the initial value of r, to zero. Choosing 
C = — (z + (1 + 2p)A d D)/(\ + p) maintains the u equa- 
tion at a fixed point with u* = u(l = 0). Then one need 
only look for stable fixed points of the D equation. 

For d > 2 (e > 0), this equation has a stable "weak- 
coupling" fixed point at D* = 0, which supplements the 
other fixed point values: r* — 0, and u* — u(l = 0). 
This fixed point gives rise to mean-field exponents (e.g., 
z = 2 and v = 1/2 follow immediately from the non- 
renormalization of the response function) , and is reached 
for initial values of D less than D c = e/A d . 

For D > D c , however, the recursion relation for D runs 
off to D = 00, rendering the critical behavior incalculable 
by perturbative techniques. Presumably the transition is 
controlled by a nonperturbative "strong-coupling" fixed 
point in this regime. For d < 2, however, the weak- 
coupling fixed point with D* = is unstable for any 
positive D, and strong-coupling behavior always obtains. 

While the strong-coupling transition cannot yet be de- 
scribed in full detail, certain aspects of it have been 
elucidated in previous work. In particular, Becker and 
Kramer g] have derived results for the amount, r c , by 
which the critical value of r is shifted away from r = at 
this transition. They find that r c < (in the Ito repre- 
sentation) for d=l and d=2. There is also a rather com- 
plete solution jl],U of the single-variable (d — 0) prob- 
lem for N = 1, where r c is also shown to be negative. 
These results allow one to make the striking prediction 
that for N = 1 the susceptibility of model ([!]) at zero 
frequency (u> — 0) and wavevector (k = 0) diverges over 
some nonzero range of values of r. To understand this, 
first recall that, as discussed above, the response func- 
tion g(k,t) in the absorbing phase doesn't undergo any 
diagrammatic corrections due to the nonlinearity, and so 
is given by the linear result g(k,t) = 9(t)e~^ k +r '> t . It 
follows at once that the susceptibility at k = u> = 0, de- 
fined as x = Jo°° dtg(k = 0, t), is infinite for all negative 
values of r that lie in the absorbing phase. The previous 
results that r c < for d < 2 (which we extend to d > 2 
below) , therefore imply the divergence of x m the entire 
range > r > r c . For d — we demonstrate explicitly 
below that this range is controlled by a fixed line with 
a continuously varying exponent. Note that the critical 
exponents z and v associated with the response function 
take their mean-field values, 2 and 1/2 respectively, at 
the point r = where the susceptibility diverges. This is 
consistent with the recursion relations for r and p in (g). 
Keep in mind, however, that the strong-coupling transi- 
tion into the active phase occurs at r(l = 0) = r c < 0, 
and so is represented by an inaccessible fixed point of the 
r recursion relation with r* = —00. 

Thus for N = 1 we are led to the schematic phase 
diagrams shown in fig. 1. In constructing the diagram 
for d > 2, we have generalized ref. Q's calculation of 
r c to dimensions larger than 2. That calculation is based 
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on mapping the computation of r c onto the quantum me- 
chanical problem of finding the lowest bound state energy 
of a potential given by the spatial correlation function of 
the noise. The analysis can easily be extended to di- 
mensions d > 2, where it is well known that the depth 
of a potential well has to exceed a d-dependent critical 
value, D c , in order for there to be a bound state. The 
existence of a bound state for D > D c implies that the 
critical value of r is shifted to r c < 0, corresponding to 
the strong coupling fixed point. When D < D c , there 
is no bound state, so r c = 0, corresponding to the weak 
coupling fixed point. These results, though obtained in 
a different way, are fully consistent with our earlier RG 
analysis, and are summarized in fig. lb. 




FIG. 1. Schematic phase diagram for model (1) with N = 1 
for (a) d < 2; (b) d > 2. Weak (strong) coupling transiti on 
occurs below (above) the multicritical point P in (b); trans 
ition is always strong-coupling in (a). Susceptibility diverges 
in the absorbing phases when r < 0, and possibly also in 
portions of the a ctive phases. 

The multicritical point P in fig. lb is controlled by the 
unstable fixed point r* , u* = u(l = 0), D* = tj Ad of Eq. 
(|), with C = ~(z + (1 + 2p)e)/(l + p). The simple scal- 
ing relation, M(r) = b^Mfrb 1 ^), for the order parameter 
M =< n > yields M ~ r@, where (3 — —C,v. Since z = 2 
and v = 1/2, we infer that = (2+ (1 + 2p)e)/2(l + p) 
for this multicritical point. This result is exact for p = 0, 
but has corrections of 0(e 2 ) for p = 1. 

We now turn to the phase diagrams of the O(N) sym- 
metric models with N > 2 and p = 1. The continu- 
ous symmetry of these models allows for two distinct ac- 
tive phases: one which preserves this symmetry and one 
which breaks it. The first of these has M =< n >= 
and Q =< q >=< (n) 2 /N >^ 0, while the second has 
both M ^ Q and Q ^ @. Since for d < 2, it is 
extremely difficult to break a continuous symmetry in a 
noisy system |HJ , we anticipate that the only active state 
is the symmetric one in this case. For d > 2 both active 
phases can occur in the phase diagram. The absorbing 
phase always occurs for sufficiently large r. 



These features can be investigated explicitly in the ex- 
actly solvable limit of N = oo. In this limit, the diagrams 
contributing to the perturbation expansion for all quan- 
tities of interest are small in number and can be easily 
summed, with the following results: 




FIG. 2. Schematic phase diagram for model (1) with 
N = oo for (a) d < 2; (b) d > 2. Symmetric (bro- 
ken-symmetric) to absorbing transition is strong (weak) cou- 
pling in (b); transition is always strong coupling in (a). 

For d < 2 there is, as anticipated, no phase that breaks 
the O(N) symmetry, i.e., for which M / 0. There is, 
however, a symmetric active phase with Q > 0. The 
value of Q is determined by DId(r + uQ) = 1 (where Id 
was defined earlier). The phase boundary r c (D) between 
the absorbing and active phases is thus determined by 
DI d (r c ) = 1, so that r c (D) - D~ 2 / e as D -> 0. Note that 
r c is positive here. Hence the susceptibility at k = uj = 
remains finite at the transition, reflecting the fact that 
< n >= in the active phase. Critical exponents are 
of course also readily calculated. For example, the order 
parameter Q vanishes like (r c — r)" with = 1. 

For d > 2 the phase diagram is somewhat more com- 
plicated, containing a broken-symmetric phase for small 
D, a symmetric active phase for larger D (see fig. 2b), 
and a multicritical point where these active phases and 
the absorbing phase meet. Again, all phase boundaries 
and critical exponents can be exactly computed. For ex- 
ample, the transition from the absorbing phase to the 
broken-symmetry active phase is mean-field-like, and oc- 
curs at the unshifted critical value r c = 0. The exponent 
governing the decay of Q at the symmetric-to-absorbing 
transition continues to assume the value = 1, the phase 
boundary being determined by the equation DId{r c ) = 1, 
which has a solution only for D > D c , with D c = 1/7^(0). 
The second-order phase boundary between the two active 
phases occurs at D = D c for all r < 0, M vanishing lin- 
early with D c — D as this boundary is approached from 
the broken-symmetric phase. The multicritical point S 
occurs at r = 0, D = D c (fig. 2b). 

While many details of these phase diagrams are doubt- 
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less special to N — oo, certain qualitative features should 
continue to apply to physically-relevant values of N such 
as N = 2. The existence of both types of active phases for 
d > 2, but only the symmetric active phase for d < 2, is, 
e.g., a general feature. For d > 2, the broken-symmetric 
active phase will continue to occur for small D, where 
the exponents for the transition to the absorbing phase 
continue to assume mean-field values, and r c = 0. The 
exponents for the absorbing-to-symmetric active phase 
transition are presumably controlled by a strong-coupling 
fixed point, and so are difficult to calculate. 

We turn next to a scaling characterization of the 
strong-coupling transition for N = 1. This is very sim- 
ilar to standard scaling theories of equilibrium critical 
phenomena, but since one of phases is absorbing, an 
additional independent exponent is required for a com- 
plete description of the phase transition. To understand 
this, consider the RG analysis of the steady-state two- 
point correlation function C(x,t) =< n(x,t)n(0,0) >. 
The rescalings of n, x and t given above yield the fa- 
miliar RG equation C(x,t,Sr) = b 2C -C{x lb,tjb z ,Srb x ' v ), 
where Sr — r c — r. First consider the equal-time case, 
t = 0. Here we obtain C(x) ~ Sr~ 2C,v 'c(x / '£) , where 
c(y) is a scaling function, and the correlation length £ di- 
verges like £ ~ Sr~ v . In typical critical phenomena, one 
can go to the critical point, Sr = 0, assume that C(x) 
approaches a nonzero value in this limit, and conclude 
that C(x) ~ x 2 ^ . Here, however, because C is identi- 
cally zero in the absorbing phase, it vanishes as Sr — > 
from the active phase. Hence in the limit Sr — > 0, C(x) 
takes the form 5r A x~( d ~ 2+n \ where A is a new, appar- 
ently independent, critical exponent, and r\ is the stan- 
dard correlation function exponent. Consistency between 
this result and the previous scaling expression requires 
c(y) ~ y-( d - 2+I '' for y « 1, i.e., A = -i^(2C + - 2 + 77) . 
This scaling law can be written more familiarly in terms 
of the order parameter exponent f3, shown above to sat- 
isfy f3 = which yields (3 = (A+i/(d-2+rj))/2. This 
is a generalization of (3 = v(d — 2 + r\)j2 1 which holds for 
ordinary equilibrium critical phenomena, where A = 0. 
In the directed percolation problem where the tran- 
sition is also into an absorbing phase, a simple graphical 
argument (which fails for the current models), shows that 
A = (3, whereupon one obtains (3 = v(d — 2 + 77) Jl2]| . 

Similar considerations allow one easily to generalize 
other scaling laws to account for the new exponent A. 
For example, the autocorrelation function C(x — 0, t) is 
readily shown to decay like Sr A t~^ d ~ 2+,1 ^ z as Sr — > 0. 
The exponent 7 governing the critical singularity of the 
equal-time correlation function (i.e., the static structure 
factor), at k — is given by 7 = v(2 — rj) — A. Note that 
our earlier results, and the scaling law relating A to j3, 
yield the following results for A at the mean-field transi- 
tion and multicritical point, respectively, for small D and 
d > 2: A = 2/(1 + p) and A = (4 + (1 + 3/j)e/2(l + p). 



Of course both these transitions have 77 = and v = 1/2. 

Finally, we verify some of the general results derived 
above in the solvable case of d — 0, N — 1. Let us first 
check the predicted divergence of the uniform suscepti- 
bility over the range < r < r c , by solving the d = 
problem in the presence of a uniform field, h, where the 
equation takes the form dn/dt — —rn — un 2+p + h + nr)(t). 
As in the case h = 0, the Fokker-Planck equation || for 
the steady-state probability distribution function P(n) 
can be solved explicitly, with the (Ito representation) re- 
suit: P(n) = ± /"dnn-ar/D-ag^fc/one-W-VDCp-i). 

here Z is a constant chosen to normalize J °° dnP(n) to 
unity. Given this expression, one readily calculates < n > 
as a function of h, thereby deriving a formula for the uni- 
form susceptibility x(/i) = d < n > jdh. In the limit of 
small h one finds that x(/i) approaches a finite limit as 
h — > 0, provided \2r/D + 1| > 1. In the band of r val- 
ues defined by \2r/D + 1 < 1, however, x{h = 0) is 
infinite, diverging like /il 2r /- D+1 l~ 1 as h — > 0. Since the 
transition from the absorbing to the active phase occurs 
jl],U & t r c = —D/2, this is consistent with our general 
argument that x(h — 0) is infinite throughout the range 
< r < r c whenever r c < 0. The exact solution shows 
that for d — the divergence actually extends into the 
active phase as well. The explicit dependence on r of 
the exponent governing this divergence as h — > implies 
that the interval \2r/D + 1| < 1 is controlled by a fixed 
line of the RG, with a continuously- varying exponent. It 
is easily shown that x{h) diverges logarithmically at the 
ends of this interval, and has logarithmic corrections at 
the critical point r = r c . 

Next we use the d = results to check the scaling 
theory developed above. In reference || one finds that 
the autocorrelation function C{t) decays like Srt~ x l 2 for 
t << t, where the characteristic time r diverges like 
Sr~ 2 for small Sr. The order parameter decays like 
Sr 1 . According to our earlier definitions, these results 
imply exponent values A = l,/3 = 1, vz = 2, and 
(d — 2 + rf)/z = 1/2, respectively. It is a trivial matter 
to verify that these values satisfy our proposed scaling 
relation 2(3 = A + v(d - 2 + 77). 
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